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Quantum entanglement, after playing a signifi- 
cant role in the development of the foundations of 
quantum mechanics [Q-Q], has been recently redis- 
covered as a new physical resource with potential 
commercial applications such as, for example, quan- 
tum cryptography [Q, better frequency standards |Q 
or quantum-enhanced positioning and clock synchro- 
nization On the mathematical side the studies 
of entanglement have revealed very interesting con- 
nections with the theory of positive maps P|J8|. The 
capacity to generate entangled states is one of the 
basic requirements for building quantum computers. 
Hence, efficient experimental methods for detection, 
verification and estimation of quantum entanglement 
are of great practical importance. Here, we propose 
an experimentally viable, direct detection of quantum 
entanglement which is efficient and does not require 
any a priori knowledge about the quantum state. In 
a particular case of two entangled qubits it provides 
an estimation of the amount of entanglement. We 
view this method as a new form of quantum compu- 
tation, namely, as a decision problem with quantum 
data structure. 

Suppose we are given n pairs of particles, all in the 
same quantum state described by some density opera- 
tor g, which is unknown. We need to decide whether the 
particles in each pair are entangled or not. From a math- 
ematical point of view we need to assert whether g can 
be written as a convex sum of product states M , 
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with | (Xif and | /Jj) pertaining to different particles in the 
pair, and ^ Pi = 1- It is assumed that the Hilbert spaces 
associated with each particle arc of finite dimensions d 
(taken to be the same for the two particles), so that one 
can always find k < d 2 . If g were known then we could 
try either to find the decomposition directly or to 
use one of the mathematical separability criteria M . For 
sufficiently large n we may indeed start with the quan- 
tum state estimation, however, this involves estimating 
d 4 — I real parameters of g, most of which are irrelevant 
in the context of the entanglement detection. In the fol- 
lowing we describe a direct method of detecting quantum 
entanglement without invoking the state estimation. 



We construct a measurement which can be performed 
on all copies of g and which is as powerful in detect- 
ing quantum entanglement as the best mathematical test 
based on positive maps [Q. The measurement can be 
viewed as two consecutive physical operations: firstly, 
we construct a transformation which maps g into an ap- 
propriate state g' and, secondly, we measure the lowest 
eigenvalue of g' . This eigenvalue alone serves as a sepa- 
rability indicator. 

A convenient starting point for our construction is the 
most powerful, albeit purely mathematical and not di- 
rectly implementable, separability criterion proposed to 
date. It is based on mathematical properties of linear 
positive maps acting on matrices Q| . Let Md be a space 
of matrices of dimension d; recall that A : Md i— > Md is 
called positive if X > implies A(X) > (expression 
X > means that the matrix X has a nonnegative spec- 
trum) . If the induced map I <£> A is also positive then A 
is called completely positive, and, as such, it represents 
a physically allowed transformation of density operators 
(here I denotes the identity map on an auxiliary system 
of any dimension). Using this terminology the separabil- 
ity criterion reads : g is separable iff 



[I<g) A] (g) > 0, 



(2) 



for all positive but not completely positive maps A : 
Md i— ► Md acting on the second particle. In fact it is 
sufficient to consider only positive maps A such that the 
maximum of Tr A(g) over all g is equal to unity. Other 
positive maps differ only by a positive multiplicative fac- 
tor which does not affect the condition (||). 

Furthermore, in some cases, instead of scanning all 
positive maps, we can choose just one. For example, all 
entangled states of two qubits can be detected by choos- 
ing A to be transposition |f(],[7j . The snag is that positive 
maps A, such as an anti-unitary transposition, and the 
induced maps I <£> A cannot be implemented in a labora- 
tory. Thus, the criterion (Q) tacitly assumes prior knowl- 
edge of g. However, there is a way to modify it so that 
it becomes experimentally viable without involving any 
state estimation. 

If we mix in an appropriate proportion [I ® A] with 
a depolarizing map that turns any density matrix into 
a maximally mixed state then the resulting map can be 
completely positive This is because the lowest negative 
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eigenvalues generated by the induced map [(I®I)®(I®A)] 
can be offset by the positive eigenvalues of the maxi- 
mally mixed state generated by the depolarizing map. 
The most negative eigenvalue — A < is obtained when 
[(I <£> I) <B> (I ® A)] acts on the maximally entangled state 

of the form -7= X^=i l*)K)> where each state | i) pertains 
to a d 2 dimensional subsystem which itself is composed 
of two d dimensional parts. Thus the map 



pr®A](<?)=p 



i®i 
~d 2 ~ 



(l-p)[l®A](g), 



(3) 



is completely positive and therefore physically imple- 
mentable when the induced map [(I ® I) ® (I ® A)] is 
positive, which happens for p > {d i X)/(d i X + 1) 



By inserting the threshold value p ■ 
into (JU) we can modify the criterion 
separable iff for all positive maps A, 



[I® A]( e ) > 



d 2 X 
d 4 X + l' 



(d 4 X)/(d 4 X + 1) 
as follows: g is 
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i.e. when the minimal eigenvalue of the transformed state 
g' = [I® A](g) is greater than (d 2 X)/(d i X + 1). In gen- 
eral, for some maps A, the related completely positive 
maps I ® A are not trace-preserving and require posts- 
elections in their physical implementations. Maps such 
as I ® A have been referred to as "structural" physical 
approximations of unphysical maps I « 



A |11 



For example, if we take A to be transposition T, (the 
first positive map used for detecting entanglement), we 
obtain 



■®T](g) = 



d 3 + 1 



I® I 



1 
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T](g). (5) 



In the two qubit case, where the partial transposition is 
a sharp test for entanglement, we obtain, 



-i ® i ■ 



T](g), 



9 9' 
which can be represented and implemented as 



1 2 

-Ai <g> A 2 + -I ® (J x o z k\u z u m , 



(6) 



(7) 



with the two channels defined as: 
A i(£0 = o ^2 <?iQVt, A 2 (e) = j a ie a 'i- ( 8 ) 



z—o,x : y 7 z 



The map can be implemented by applying selected 
products of unitary (Pauli) transformations with the pre- 
scribed probabilities. The map It is trace-preserving 
hence any postselection in experimental realizations is 
avoided. 

Thus, in order to detect entanglement of an arbitrary 
two-qubit state g it is enough to estimate a single pa- 
rameter, i.e. the minimal eigenvalue of \l®T](g). The 



state g is separable iff this eigenvalue satisfies X m i n > | . 
Let us also point out an extra bonus: X m i n gives us —A', 
the most negative eigenvalue of [I ® T](g), which enters 
the expression for the upper and lower bounds for the 
entanglement of formation, 



< H 
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where H(x) is the Shannon entropy. The above formulae 
can be derived from the estimations of the concurrence 
provided by Verstraete et al |l2| . 

Suppose for a moment that I ® A is trace-preserving, 
e.g. the transposition case. The first part of our en- 
tanglement detection measurement is accomplished by 
applying I ® A to each of the n pairs to obtain n copies 
of g' = [I® A] (g). Then, following the criterion (^), we 
need to measure the lowest eigenvalue of g' . 

This can be viewed as a special case of the spectrum 
estimation and possible approaches depend a lot on par- 
ticular physical realizations of g' . Here, we provide two 
general solutions. The first one, based on quantum intcr- 
fcrometry, is conceptually simple and relies on estimating 
d 2 — 1 parameters from which the spectrum of g' can be 
calculated (this is a significant gain over the state esti- 
mation which involves d 4 — 1 parameters). The second 
solution is a joint measurement on all copies of g' which 
gives directly the estimate of the lowest eigenvalue. 

We start with the quantum interferometry, presented 
here as a quantum network shown in Fig.([l]). A typical 
interferometric set-up for a single qubit — the Hadamard 
gate, phase shift <j>, the Hadamard gate, followed by a 
measurement is modified by inserting in between the 
Hadamard gates a controlled-^ operation, with its con- 
trol on the qubit and with U acting on a quantum system 
described by some unknown density operator p. (N.B. 
we do not assume anything about the form of p, it can, 
for example, describe several entangled or separable sub- 
systems.) The action of the controlled-C/ on p modifies 
the interference pattern by the factor, 
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where v is the new visibility and a is the shift of the 
interference fringes, also known as the Pancharatnam 
phase [Q. Formula ( |l0| ) has been derived, in the context 
of geometric phases, by Sjoqvist et al. Q. 
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FIG. 1. Both the visibility and the shift of the interfer- 
ence patterns of a single qubit (top line) are affected by the 
controlled-17 operation. This set-up allows to estimate Tr Up, 
the average value of U in state p. 

The network can evaluate certain non-linear function- 
als of density operators. Indeed, let us choose p to be 
composed of two subsystems, p = g a ®gb, and let U to be 
the exchange operator V such that V a) 0) — \0)\ ot) 
for any pure states of the two subsystems. The inter- 
ference pattern is now modified by the factor Tr V(g a ® 
Qb) = Tr Q a Qb- For p = g ® g we can estimate Trg 2 , 
which gives us an estimate of $^i=i ^? i where Xi are the 
eigenvalues of g. N.B. Tr g 2 is real hence there is no need 
to sweep the phase if in the interferometer, it can be fixed 
at <p = 0. 

In general, in order to calculate the spectrum of any 
m x m density matrix g we need to estimate m — 1 
parameters Trg 2 , Trg 3 ,... Tig m . For this we need the 
controllcd-shift operation. Given k systems of dimension 
m we define the shift as 

vW\4>i)\4>2)-\4>k) = \<t>k)\cl>i)-\<t>k-i), (ii) 

for any states | 0) . Such an operation can be easily con- 
structed by cascading k — 1 swaps V. This time, if we 
prepare p = g® k the interference will be modified by the 
factor 
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Thus measuring the average values of for k = 

2,3...m gives us effectively the spectrum of g. In par- 
ticular, in our case, we obtain the spectrum (and the 



lowest eigenvalue) of g' = [I ® A](f>) by estimating d 2 — 1 
parameters Trg' k , where k = 2...d 2 . Again, the phase in 
the interferometry can be fixed at tp = 0. 

The interferometric scheme described above is concep- 
tually simple and experimentally viable, however, if the 
simplicity of the implementation is not an issue then we 
can measure the estimate of the lowest eigenvalue di- 
rectly. This requires a join measurement on all of the n 
pairs. We use the Keyl and Werner spectrum estimation 
method [jl5|| , which, in the entanglement detection con- 
text, works as follows. The n copies of the mxm state g' 
(in our case m = d 2 ) form an operator on the n-fold ten- 
sor product space which can be decomposed according to 
irreps of SU(m), so that each summand, including mul- 
tiplicities, is labelled by a Young tableau, i.e. n boxes 
arranged in rows of decreasing length (c.f. [[l6| for the 
SU{2) case). The tableaus give a family of projectors for 
the spectrum estimation measurement. The normalized 
row lengths of each tableau are taken as estimates of the 
ordered sequence of eigenvalues of g' . The probability 
that the error is greater than some fixed e decreases ex- 
ponentially with n [ p^[ . In our particular case, we are 
interested only in the lowest eigenvalue. We modify the 
Keyl- Werner scheme by adding together all projectors 
corresponding to Young tableaus with the fixed length 
of the last row. The measurement determined by these 
projectors gives directly the estimate of only one param- 
eter — the lowest eigenvalue of g 1 . Such a measurement 
can be represented as a quantum network implementing 
projections on the symmetric and on partially symmet- 
ric subspaces (see JTtJ for the network projecting on the 
symmetric subspace). 

Our considerations remain valid, with some minor 
modifications, when I (g> A is not trace-preserving. In 
this case experimental implementations require postse- 
lections, which result in n' = nTr (I £g> A(g) copies of nor- 
malized states I £g> A(g)/Tr (I £g) A(g)). The spectrum es- 
timation procedure is not affected, however, before check- 
ing the condition (Q) the lowest eigenvalue has to be 

rescaled by the factor Tr (I £g> A(g)). 

Let us summarize our findings. Given n copies of a bi- 
partite d®d system described by some unknown density 
operator g we can test for entanglement either by esti- 
mating g and applying criterion (|2|), or, more directly, 
by performing the measurements we have just described. 
The state estimation involves estimating d A — 1 param- 
eters of g, most of which are of no relevance for the en- 
tanglement detection. The optimal state estimations rely 
on joint measurements on all copies of g. However, one 
can construct less efficient but simpler state estimation 
methods which involve measurements only on individual 
copies. Our more direct, interferometry based, method 
requires estimations of only d 2 — 1 parameters and joint 
operations on d copies of g 1 . The most demanding, from 
the experimental point of view, is our second method. It 
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is a measurement with an outcome which is an estimate 
of just one parameter, but, like the optimal state esti- 
mation, the measurement involves joint operations on all 
copies of g' . Both direct and indirect entanglement detec- 
tions have their own merits. Depending on the context, 
applications, and technologies involved one can choose 
one or the other. 

Direct entanglement detections, can be employed as 
sub-routines in quantum computation. For example, one 
may consider performing or not performing a quantum 
operation on a given quantum system conditioned on 
some part of quantum data being entangled or not. In 
fact direct entanglement detections can be viewed as 
quantum computations solving an inherently quantum 
decision problem: given as an input n copies of g de- 
cide whether g is entangled. Here the input data is 
quantum and such a decision problem cannot even be 
even formulated for classical computers. Nonetheless the 
problem is perfectly well defined for quantum comput- 
ers. Finally, let us add that the method presented here 
can be easily generalized to cover all linear maps tests 
for arbitrary multiparticle entanglement Jl8[ ] and the so 
called k-positive map tests detecting Schmidt numbers 
of density matrices Jl9| . Modification of the method to 
the case of two distant labs (i. e. under restrictions to 
local operations and classical communication ||]) will be 
considered elsewhere. 

For the sake of completeness we should also mention 
here that there are two-particle observables, called en- 
tanglement witnesses which can detect quantum entan- 
glement is some special cases (see [^0]||). They have 
positive mean values on all separable states and nega- 
tive on some entangled states. Therefore any individual 
entanglement witness leaves many entangled states un- 
detected. When g is unknown we need to check infinitely 
many witnesses, which effectively reduces this approach 
to the quantum state estimation. However, let us point 
out any witness defines a positive map which can be used 
in our test. 

To conclude, we have demonstrated that direct and 
physically implementable methods of entanglement de- 
tection are possible. They are equivalent to the most 
powerful mathematical separability criteria known to 
date M. We have described their possible realizations 
in the generic terms of quantum gates and networks. 
These generic components can be implemented using sev- 
eral experimental techniques ranging from trapped ions 
to quantum dots pT| . 
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